DISTANCE fc-GRAPHS OF HYPERCUBE AND g-HERMITE 

POLYNOMIALS 



HUN HEE LEE AND NOBUAKI OBATA 

Abstract. We will prove that some weighted graphs on the distance fc-graph 
of hypercubes approximate the g-Hermite polynomial of a g-gaussian variable 
J^jfy by providing an appropriate matrix model. 



1. Introduction 

During the last decade quantum probabilistic approach to the study of asymp- 
totic spectral analysis of graphs has been developed considerably, and many impor- 
tant probability distributions are obtained as scaled limits of spectral distributions 
of growing graphs, see [7J and references cited therein. For example, the spectral 
distribution of the Hamming graph H(d,n) converges to the standard normal dis- 
tribution or the Poisson distribution by taking a suitable scaling limit as d, n — > 00, 
which was first proved by means of quantum probabilistic method in [6]. 

In this paper we extend quantum probabilistic method to weighted graphs and 
derive the distribution of the g-gaussian variable G q (—1 < q < 1) by Bozejko and 
Speicher [4]. In fact, we will show that the probability distribution o{H%(G q ), where 
H% is the fc-th g-Hermite polynomial, is derived from a sequence of weighted graphs 
on the distance fc-graph of n-cube. It is noteworthy that g-gaussian variables for 
q ^ 0,1 have not been so far observed in line with the asymptotic spectral analysis 
[TJ, where consideration has been mostly restricted to simple graphs, i.e., undirected 
graphs with no loops and no weights (multi-edges). While, our result in a particular 
case of q = 1 is proved in [9] in an alternative manner. 

This paper is organized as follows. In section[2]we will recall some basic notations 
and tools we need for the proof of the result, which includes the baby Fock model 
for g-gaussians. In section [3] we will explain the connection between the distance 
fc-graph of hypercubes and the baby Fock model, and finally we will prove the main 
result in section |4j For basic concepts of quantum (non-commutative) probability, 
see e.g., [7j[8]. 



2. Preliminaries 

2.1. Non-commutative probability and convergence of distributions. Let 

(A4 n , p n )n>i and (M.,<p) be non-commutative (tracial) PF*-probability spaces. For 
to > 1 we consider self-adjoint elements X", • • • , g A4 n and X\, ■ ■ ■ , X m E A4. 
We say that the sequence of n-tuples {X™, • • • , X^) converges to (X±, ■ ■ ■ , X m ) in 
distribution as n — >■ 00 if 

^(P(A7,-.. ,XI))^^{P(X X ,--- ,X m )) 

for any polynomial P of to non-commuting variables. We will use the notation 

/ vn \rn \ dist / v v \ 
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Remark 2.1. In classical probability theory, we say that a sequence of random 
variables X n converges to X in distribution if the probability distribution of X n 
converges weakly to that of X . Therefore the above notion of convergence in distri- 
bution for non-commuting random variables is not compatible to the one in classical 
probability theory. We say often that a sequence of random variables X n converges 
to X in moments if the m-th moment of X n converges to that of X for all m. It is 
known that the convergence in moments implies the convergence in distribution if 
the distribution of X allows a positive answer to the determinate moment problem. 
For more information along this line, see e.g., [5) Chapter 4]. 

We present a lemma, which is probably a folklore in non-commutative probability. 
We include the proof for the convenience of the readers. Recall that the non- 
commutative L p -space L P (M, ip) is the completion of (Ai, \\-\\ p ), where the p-norm 
on M 1 H'll is defined by 

\\A\ P = W))', 

Lemma 2.2. Let X", • • • , X™ 6 M n and X\,--- ,X m 6 M be self-adjoint ele- 
ments and assume that 

( vn x^n \ dist , v v \ 

{^1 ) " ' ! A mJ > 1^1) " " ' i-^ra)- 

Let Yp, ■ ■ ■ ,Y£ € M n be self-adjoint elements satisfying 

(1) {||X"|| , ||^i"|L '■ n > 1,1 < i < m} are bounded for all p > 1 with bounds 
depending only on p. 

(2) \\X? - Y^\\ p -^Oosn^oo for all p> 1. 

Then, we also have 

(YT,-- - ,Y£)2$(Xi,.~,X m ). 

Proof. It is enough to prove the case of monomials. Let us fix a monomial P. Our 
strategy is the following. We begin with lim„_ J . 00 ip. n (P{X™ , • • • , X™ )) and exchange 
each X™ with Y™ one by one. Here we use the Holder inequality and the conditions 
(1) and (2). Then we replace X^ by Y£ similarly, and we repeat the same procedure 
until we replace all X™'s by Y™'s. 

Let us do the case of (p((X™) 2 (X^) 3 ) as an example. Then, by the Holder 
inequality we have 

\^(x?-Ynxux2) 3 )\ < \\x?-Y?\\jx?\\jxnh 

so that the conditions (1) and (2) allow us to conclude that 

lim V ((X?f{X%Y) = lim ^X-{X-f). 

n—^oo n—^oo 

Note that the number 5 in ||-|| 5 is the degree of the monomial. 

□ 

2.2. The g-gaussians. In this subsection we recall the q-gaussian variables. Let 
JTr be a real Hilbert space and H = TJr + iHs. be its complexification. We consider 
the operator of symmetrization P n on %® n defined by 

p n = si, p n (h ® •• • ® /„) = ® • •• ® Aw, 

where S n denotes the symmetric group of permutations of n elements and 

= #{(*:J')I 1 < hj < n,ir(i) > vr(j)} 

is the number of inversions of tt € S n . 

Now we define the q-inner product (•, ■} on T-i® n by 

(t,v) a = {t,PnTl) for tveU® n , 



where (•, •) is the inner product in *H® n . When —1 < q < 1 P n 's are strictly positive 
so that (*; m ) q i s actually an inner product. We denote by 'U^i 71 the resulting 
Hilbert space. Then one can associate the q-Fock space F q (H). 

n>l 

where f2 is a unit vector called vacuum. When q = we recover the so-called full 
Fock space over T-L. In the extreme cases q — ±1, we have 

n>l n>l 

which are referred to as the Bosonic and Fermionic Fock spaces, respectively. Here 
® s and ® a refer to symmetric and anti-symmetric tensor product of Hilbert spaces, 
respectively. 

For h € .Hr, we can define a generalized g-semi-circular random variable on 
F q {%) (-K q < 1) by 

s(h)=£ q (h)+t q (h), 

where £ q (h) is the left creation operator by h 6 H and ^(/i) is the adjoint of £ q (h). 
Recall that £ q (h) is given by 

£ q (h)fl = h, £ q {h)fx ®---®f„ = h®f l ®---®f n . 

Note that in the extreme cases q = ±1 the creation operators have slightly different 
forms due to the symmetrization procedure. For example, we have 

1 - 

£x(h)Q = h, £Ah)f® n = V f® n ~ k ®h® f m . 

n + 1 ^-^ 

Then, we get a -probability space (T q , ip q ), where T q is the von Neumann algebra 
generated by {s(h) : h £ Hr} in B(J r g ('H)) and <p q (-) = { ■ £1, 0) is the vacuum state. 

2.3. Baby Fock space, central limit procedure and hypercontractivity. In 

this section we collect background materials focusing on baby Fock space and its 
central limit procedure due to Biane [2]. 

Let / = {1, 2, • • • , n} be a fixed index set and e : I x / — > {±1} be a "choice of 
sign" function satisfying 

e («, j) = e(M) = -1) e /• 

Now we consider the unital algebra A{I,e) with generators (xi)i^i satisfying 

In particular, we have xf = 1, i & I, where 1 refers to the unit of the algebra. The 
involution on A(I,e) is defined by x* — Xj. We will use the following notations for 
the elements in A(I, e). 

x := I and xa ■= x h ■ ■ ■ x ik , A = {h < ■ ■ ■ < i k } C I. 

Then, {xa ■ A C /} is a basis for A(l,e). Let : A(I,e) — > C be the tracial state 
given by 

4> e {x A ) = S A< 0, 

which gives rise to a natural inner product on A{I, e) as follows: 

(x,y) := (f> £ (y*x), x,yeA(I,e). 

Let 

H = L 2 (A(I,s),^) 

be the corresponding L 2 -space, then clearly {xa '■ A C 1} is an orthonormal basis 
for ff. 
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Now we consider left creations P* and left annihilations fa in B(H), i £ / in this 
context. 

fXiXA if i i A XiXA Hi £ A 

, Pi(x A )={ , i £ I, AC. I. 

if i £ A [ if i£ A 

The baby gaussians are defined as follows. 

7l := + P u l<i<n. 

It is straightforward to check that ji is the same as left multiplication operator with 
respect to Xi , so that we often identify 7, and Xi . Let T n be the von Neumann algebra 
generated by {7i}"=i i n B{H). It is also known that 1 is a cyclic and separating 
vector for r n , and the above is the faithful GNS representation of (T„,t^), where 
is the vacuum state on T n given by 

<(•) = <• 1,1). 

Note that is tracial on T n . 

Now we replace {e(i, j) : 1 < i < j < n} with a family of i.i.d. random variables 
with 

P( £ (*,. ? ) = -l) = I_i, P(e(i,j) = 1) = I±*. 

We set 

1 - 

Then, the Speicher's central limit procedure ((TU] or [5]) tells us the following. 
Proposition 2.3. We have 

dist ^ 
S n r Kjj q 

for almost every e. In other words, we have for any polynomial Q 

hm T*(Q(s n ))=T q (Q(G q )) 

n—too 

for almost every e. 

We close this section with the last ingredient for the proof, namely hypercontrac- 
tivity of baby Ornstein-Uhlenbeck semigroup. Recall that the number operator N e 
on H is given by N e = J2i^i fit Pi- Then, for any A C I we have N e XA = \A\ xa- 
Since 1 € H is separating and cyclic for T n we define the e-Ornstein-Uhlenbeck 
semigroup F t e : T n -» T n by 

(2.1) P t £ (X)l = e- tN '(Xl), X£T n . 

Theorem 2.4. Let 1 < p < r < 00 we have 

- J»-l 



\ p t\\Lv^LT < 1 */ and onl V if e " 
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The above result is due to Biane [5J Theorem 5], which leads us to a Khinchine 
type inequality as follows. 

n 

Corollary 2.5. Let X = 

a ii—ifcTii ' ' ' 7ifc £ then we have 

il,— ,i fc = l 

||X|| 2 <||X|| p <(p-l)t||X|| 2 . 
Proof. Note that P t £ (X) = e _fct X. Then it is a direct application of Theorem 

nn □ 



3. The distance /c-graph of hypercubes and Baby Fock model 



We begin this section by recalling some graph theoretic notions. A weighted 
graph consists of the set of vertices V, the of edges E C V x V and a collections 
of weights A = {dxy)(x,y)eE, a xy € K. We say that (V,E,A) is a weighted graph 
on the unweighted graph (V,E) (i.e. every weight is identically 1). The matrix 
A is called the adjacency matrix of the weighted graph. When the graph is finite 
(i.e. finite number of vertices) A is a non-commutative random variable in a Im- 
probability space (MiyhV'e)j where M\ v \ is the matrix algebra of the size |V| and 
(fie(B) = (B8 e , S e ) = B ee for a fixed point e <G V and B G Miyi. 

The n-cube is the unweighted graph KV, — K 2 x • • • x K 2 , the n-fold direct product 
of the complete graph K 2 with two vertices. Then there is a bijection between the 
vertices of K 2 and the set of pre-described orthonormal basis of H as follows. 

where <E {0, 1} and x\ — 1. Note that x^ 1 ■ ■ ■ x r ™ — xa with A = {1 < i < n : 

n = l}. 

In n-cube there is an edge between two vertices (?r, • • • , r n ) and (si, • • • , s n ) if 
and only if there is only one index i such that r» ^ Sj. T/ie distance k-graph of K 2 
(1 < fc < n) is the graph with the same set of vertices but with the set of edges 
described as follows. 

{((n, • • • , r n ), (si, ■ ■ ■ , s n )) : there are exactly k indices 1 < i < n with 7^ Sj}. 

Now we focus on the the Bosonic case for the moment, which implies that q = 1 
and so that the choice of sign function e is identically 1. Then, it is straightforward 
to see that the operator 

n 
i=l 

represents the adjacency matrix of K 2 since 7» is nothing but multiplying 2^. It is 
again straightforward to see that the operator 

j£[ 7»i 7i fc = 7»i 7*2 " ' 7u- 

tir" »i<»2< - "<»fc 

represents the adjacency matrix of the distance fc-graph of . Here the summation 
in the left-hand side is taken over distinct i\, ■ ■ ■ ,i^ taken from {1, • • • , n}. 

Now let us turn our attention to the case of general e. This time 7$ still acts 
on x^ 1 ■ ■ ■ x 7 ^ 1 by multiplying Xi on the left, but we need to take the commutation 
relations into account. Thus, the operator 

n 

i=i 

corresponds to the adjacency matrix of the weighted graph on K 2 described as 
follows: We put the weight 

e(r,i) = e(i,l) ri ■■■e(i,i- l)^- 1 

on the edge between r = (n, • • • , r„) and s = (si, • • • , s„) with Sj 7^ rj, but Sj = rj, 
Vj 7^ i Similarly, the operator 

E 7*1 ■■ ■ 7» fc 

represents the adjacency matrix of a weighted graph on the distance k-graph of K 2 . 
The rule for assigning weights on each edge is the following. If there is an edge 
between r = (?r,--- , r n ) and s = (s\,--- ,s n ), then we have k distinct indices 

5 



{ji < ■■■ < jk} with s jl ^ r jt , 1 < I < k but Sj = r , Vj <£ {j u ■ ■ ■ ,j k }. For a 
permutation a G S k and j = (ji, • • • , jk) we define e(cr,j) to be the number given 

by 

On the above edge we put the weight 



1 fe 



fc! 



4. Weighted graphs on Hypercubes and <?-Hermite Polynomials 

In this section we focus on the analytic part of our results, namely the convergence 
analysis. 

We recall the normalized q-Hermite polynomial H% given by the following recur- 
rence relations. 



H q (x) = 1, Hf(x) = x 

xHl(x) = Ht +1 {x) + \k] q Hl_^x), k > 1, 



where [k] q = 1 T ^ q -- 

We are interested in the following operator. 



X n ,k '■— k 7n 



■■■7ik- 

We would like to find a recurrence relation regarding X„ ; fc by multiplying X„^. 
Now we have 



J2 T*i 



■ ■ ■ 7i* 



= E 7ii ' ' ' 7 4fc 7» + E 7ii ' ' ' 7i*7* + ■■■+ E 7ii ' ' ' 7i*7i- 

*ir" *=*ir"i** *ir"i**=* 

By relabeling we denote ^ 7u ' ' ' 7i fc 7i = E 7il ' ' ' 7**+i ■ For tlic 

*i) •••>**>* ii,-- ,**i»*+i 

ond term we have 

y ] 7»i ' ' ' 7'ifc7i = ^ ] e i,i2 ' ' ' £ M*7*2 ' ' ' 7ifc 



sec- 



=*1,'" M2,'" i** 

5^ / 



^ ' £ Ml ' ' ' £ i,ifc_i7«i ' ' ' 7i*-i 



by relabeling again. If we repeat the similar relabeling, then we get 
(4.1) 




+ Y • • ■ + £ Mj • • • e M*-i + ■ • • + 1 )7i l ---7» fc _i- 

ill"' Ak-l ,-*■ ,»h-l) 

Now we define 

Y n ,k-=\ Y] [k + l^Zfa, ■ ■ ■ ,i k hh ■ ■ -H k , 

where 

:=- (£»,»! • • • £i,i h + e»,i 2 • • • £t,» fc H hi). 

Here comes the key recurrence relation obtained by dividing (|4.ip with n~^~ . 
(4.2) X n , fe+ i =X n , fc X„, 1 - [fc] g r n , fe _i, fc>l. 

We will use Lemma 12.21 for the sequences (^C n; i, • • • . Xn.k* Yn.\ : ' ' ' ^^n,k— l) and 
(X„.i, • ■ ■ , X n< k, X n< ±, • • • , X n) fe_i). In order to do so we need to check the condi- 
tions (1) and (2) in this case. 

The condition (1) can be easily checked by Corollary [53] Note that 

{lh ' • "7i* : an distinct , ik} 

is an orthonormal family in L 2 (T n ) and the number of indices (ij)j =1 with all 

distinct entries is strictly smaller than n? . 

The condition (2) is more involved. We first need to understand the limit behav- 
ior of Z(i\, • • • , ik) as n —¥ oo. Note that the sequence of random variables 

{Wi = e Ml • • • e i>ih + e iii2 ■ ■ ■ e i<ih H h 1 : i ^ ■ ■ ■ , ik-i)} 

is an independent collection. Indeed for i ^ i' both of them are not in (i±, ■ ■ ■ , ik—i), 
we have ^ (i',ij>) and ^ for any 1 < j.j' < k. Moreover, each 

random variable has mean 

nWi) = q k +q k - 1 + --- + l = [k + l] q 

and variance 

k 

a 2 = E(W?) - E{W t ) 2 = k + 1 + Y 2 (k + 1 - j)q k ~ [* + 1],- 

3=1 

Thus, the classical central limit theorem tells us that 

yfc([k + l]- x Z(ii, • • • ,i fc ) - 1) ^ AA(0, a 2 ) 

in our sense (i.e. convergence in moments, see [3] section 30] or [I]), so that their 
L p -norms also converge. Thus, we have 

\\yfc([k+l]- 1 Z(i 1 ,--- ,i fc )-l)|| p -> \\M(0,a 2 )\\ p = a p , 
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and consequently 
(4.3) 



[[k + i^zin,..- ,i k )-i\\ <^±1 

H Up n 



for big enough n. Note that a p is a constant depending only on k, q and p, so that 
it is independent of n. 
Now we get the estimate 



(4.4) 



for big enough n. Indeed, we have 



I? 2 



,«fc))7n •••7i, 



and by Corollary 12.51 Minkowski's inequality and (|4.3|) we have 

||x„, fe -r n , fc || p = [E(||x^-r n £ )fe || p 



< 



< 



< 



(p-i) 

n 2 



(p-ir 



E 



J2 llk + l^ZiH,--- ,ik)-l\ 



l I'll"' i?>k 



J2 ntfc + i]- 1 ^!,--. ,*fc)-l||* 

v * 



(p-l) = (a p + l) 



where fc is the value of X n ^ for a fixed choice of e. A standard application of 
Borel-Cantelli lemma with (14.41) leads us to the following conclusion. 

\\Xn,k ~ ^n,fe|| p 

for almost every e. This observation tells us, together with Proposition ^. 3[ that we 
can choose a specific "choice of sign" e satsfying 

dist 



Gn 



-> 0, Vp > 1 



Note that we used the fact that the collection of all monomials in non-commuting 
variables is countable. From now on we will fix this choice of e and by abuse of 
notation we will denote fe and Y£ k simply by X n ,k and Y n ,k, respectively. This 
explains how we get the condition (2). 

Finally, we present the main convergence result. 

Theorem 4.1. We have 

(X n ,l,''' i X n ,k, Yn,l' " ' i^n,fc— l) 



(Hf(G q ), ■ ■ ■ ,H%,(Gq),Hl(Gq),--- ,Hl_ x {G q )). 



Proof. We will use induction on k. When k = 1 we get the result directly from 
Proposition ^. 31 When fc = 2we note that 

2 

-2 T _ rjq. 



x n ,2 = - v 7i7j = [ —j= 

n frf V v n 



-I = X£ A -I = H«(X ntl ). 



Thus, we clearly have 

(X nA ,X rh2 ,X 7hl ) ^4 (H*(G 9 ),H*(G 9 ),H*{G q )). 

If we compare two sequences (X ni i, X Hi 2, -i^n.,1) ancl (-^n,ii -^n,2j then we get 

the result we wanted by Lemma 12.21 since we already checked the conditions (1) and 
(2) are satisfied. 

Now we suppose the the conclusion is true upto k and let us check the case for 
k + 1. Recall the recurrence relation 

X n ,k+1 = X n ^X n ^i — [k] q Y n! k-l, 

so that we get 

n (Hf(G q ),--- ,H«(G q ),H« +1 {G q ),H?(G q ),--- .^(G,)). 
This can be trivially extended to the following convergence. 

,H«(G q ),H« +1 (G q ),H?(G q ),--- ^^(G^H^G,)). 

If we compare two sequences (X„ a ,--- , F n> i, • • • ,Y n , k ^ l7 X n , k ) and 

(X ni i,--- ,X„ ife ,X n , fc+ i,y„,i, • • • ,Y n>k -i,Y n , k ), then we get the result we wanted 
by Lemma 12.21 since we already checked the conditions (1) and (2) are satisfied. 
This finishes the induction process. 

□ 

Corollary 4.2. For any k > 1 we have 

X n , k ^4 Hl{G q ). 

Thus, we can conclude that H^(G q ) can be approximated (in distribution) by a 
sequence of weighted graphs on the distance k-graph of K% . 

Remark 4.3. (1) When g = lwe recover the result of N. Obata in [9]. 

(2) The above result provides a matrix model for H^(G q ) which is homogeneous 
in degree. 
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